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^ ■ Abstract 

We discuss corrections to the Yukawa matrices of the Standard Model (SM) fcrmions in 
^ intersecting D-branc models due to four-point interactions. Recently, an intersecting D- 

O , brane model has been found where it is possible to obtain correct masses and mixings for 

all quarks as well as the tau lepton. However, the masses for the first two charged leptons 
Cn ' come close to the right values but are not quite correct. Since the electron and muon are 

\^ I quite light, it is likely that there are additional corrections to their masses which cannot be 

^vq , neglected. With this in mind, we consider contributions to the SM fermion mass matrices 

. I from four-point interactions. In an explicit model, we show that it is indeed possible to 

f^ ' obtain the SM fermion masses and mixings which are a better match to those resulting from 



experimental data extrapolated at the unification scale when these corrections are included. 
These corrections may have broader application to other models. 



'^ ■ PACS numbers: ll.lO.Kk, ll.25.Mj, ll.25.-w, 12.60.Jv 



INTRODUCTION 



In recent years, intersecting D-brane models, where the chiral fermions arise at the intersections 
between D6-branes (Type IIA) in the internal space [l| with the T-dual Type IIB description in 
terms of magnetized D-branes (^ haveprovided an exciting approach towards constructing semi- 



realistic string vacua (for reviews, see 



HQ). 



Indeed, such models provide promising setups which 
may accommodate semi-realistic features of low-energy physics. Given this, it is an interesting 
question to see how far one can get from a particular string compactification to reproducing the 
finer details of the Standard Model (SM) as a low-energy effective field theory. 

The Standard Model has an intricate structure, with three-generations of chiral fermions which 
transform as bifundamental representations of SU{3)c x SU{2)l x U{1)y- In addition to the fact 
that the SM fermions are replicated into three distinct generations, the different generations exhibit 
a distinct pattern of mass hierarchies and mixings. Interestingly, intersecting D-brane models may 
naturally generate the SM fermion mass hierarchies and mixings, as well as an explanation for the 
replication of chirality. In short, D6-branes (in Type IIA) fill four-dimensional Minkowski space- 
time and wrap 3-cycles in the compact manifold, with a stack of A^ D6-branes having a gauge group 
U{N) (or U{N/2) in the case of T^/(Z2 x Z2)) in its world volume. The 3-cycles wrapped by the 
D-branes will in general intersect multiple times in the internal space, resulting in chiral fermions 
in the bifundamental representation localized at the intersections between different stacks. The 
multiplicity of such fermions is then given by the number of times the 3-cycles intersect. 

The Yukawa couplings in intersecting D6-brane models arise from open string world-sheet in- 
stantons that connect three D6-brane intersections [5:]. For a given triplet of intersections, the 
minimal world-sheet action which contributes to the trilinear Yukawa couplings is weighted by a 
factor exp(— A(jbc)> where Aahc is the world-sheet area of the triangle bounded by the branes a, 6, 
and c. Since there are several possible triangles with different areas, mass hierarchies may inher- 
ently arise. The Yukawa couplings depend on both the D-brane positions in the internal space 
as well as on the geometry of the underlying compact manifold. Effectively, these quantities are 
parameterized by the vacuum expectation values (VEVs) of open and closed-string moduli. 

Despite substantial progress in constructing semi-realistic vacua with intersecting D-branes, 
there are many phenomenological challenges remaining, besides the usual moduli stabilization 
problem. Typically, some or all of the Yukawa couplings are typically forbidden by selection rules 
which arise from global C^(l)s which become massive via a generalized Green-Schwarz mechanism. 
It has been found that only Pati-Salam models can have all the SM fermion Yukawa couplings 



present at the stringy tree level, although some couplings which are perturbatively forbidden may 
be generated in principle via D-brane instanton effects [a, lD, Isl, Isl, lid]- Also, there has generally 
been a rank one problem in the SM fermion Yukawa matrices, preventing the generation of masses 
and mixings for the first two families of quarks and leptons. For the case of toroidal orientifold 
compactifications, this can be traced _to_ the fact that not all of the SM fermions are localized at 

lid . llJ]. However, one example of an intersecting D6- 



intersections on the same torus 



11 



12 



brane model in Type IIA on the T^/(Z2 x Z2) orientifold has recently been discovered in which 



these problems may be solved J15l . 



Id ]. Thus, this particular model may be a step forward to 



obtaining realistic phenomenology from string theory. Indeed, as we have recently shown [17], it is 
possible within the moduli space of this model to obtain the correct SM quark masses and mixings, 
the tau lepton mass, and to generate naturally small neutrino masses via the seesaw mechanism. 
In addition to these features, the model exhibits automatic gauge coupling unification, and it is 
possible to generate realistic low-energy supersymmetric particle spectra, a subset of which may 
produce the observed dark matter density. 

In spite of the successes of this model, the electron and muon masses come close to the right 
values, but they are still not quite correct. Since the electron and muon masses are very light, it 
is likely that there are additional corrections that must be considered, as first suggested a long 
time ago in [18]. This idea was applied later in string theory, for the case of the free-fermionic- 
formulation, where the general rules for calculating higher order corrections, from multipoint- 
functions, to the Yukawa couplings were first given in [19]. With this in mind, in this paper 
we consider contributions to the Yukawa couplings from four-point functions. Four-point string 
interactions in intersecting DG-brane scenarios can be regarded as the scattering amplitudes of 



four string matter fields 



201 ] ■ and are consistent to the Kahler metrics of twisted matter (211] . 



Both the quantum and classical contributions have been studied in two special cases where one 
las only an independent angle and the other has two independent angles in a closed quadrilateral 
221 ]. This was then promoted to the generalized calculation of A^-point amplitudes and four-point 



functions without angle constraints ^]. T 
closed string interactions on orbifolds 



lis analysis is based on the discussion of the twisted 



M 



261 ]. where the analogy between open strings at 



brane intersections and closed strings on orbifolds is known. The complete amplitudes for Yukawa 
couplings, as a three-point limit of the four-point amplitude from conformal field theory, are also 



discussed in 



2C 



22 



23!]. 



It is natural to ask whether the four-point function corrections to the Yukawa couplings can 
generate the correct electron and muon masses in the model discussed in 17|; indeed the correct 



electron and muon masses can be obtained by including four-point function corrections in principle. 
However, the D6-branes in this model can not form a closed quadrilateral for the required four- 
point couplings on the relevant two torus. Two stacks of D-branes overlap so the "angle" at their 
"intersection" turns out as a straight line. Instead of using the conventional calculational procedure 
in this model, we can still apply the idea of four-point function corrections by taking the vector-like 
fields at the "intersection" on the overlapped D6-branes in a asymptotic limit of the "angle" to vr, 
a physical mechanism which has not yet been clarified. Therefore, we shall consider another model 
where the standard four-point function corrections can be calculated. As we shall see later in this 
model, if the SM fermion mass matrices receive contributions from only three-point functions, one 
can only obtain the correct quark masses and tau lepton mass, but it is not possible to explain the 
CKM quark mixings, and the electron and muon masses. Introducing corrections from four-point 
functions, one can indeed obtain the correct quark and the tau lepton masses, and the CKM quark 
mixings, and as well as the electron mass. However, the muon mass is still about 36% smaller than 
the desired value. 

In this paper, we review the construction of four-point functions, and include both trilinear 
Yukawa couplings and four-point interactions to obtain the SM fermion mass matrices in an explicit 
intersecting D-brane model. We show that it is indeed possible to better match the masses and 
mixings for all fermions when four-point corrections are included. In addition, these corrections 
may have broader application to other models. 

II. GENERAL FOUR-POINT FUNCTIONS 

Let us begin by considering an open string stretched between two D-branes intersecting at an 
angle irB. From the boundary conditions we can write the mode expansion as 



dX{z) = ^ak-e- 

k 

dX{z) = Y,(^k^ez-''-'-\ (1) 



k 
where z is the worldsheet coordinate. By comparison of this expression with the mode expansion 



for a closed string in the CFT analysis of a Zjv orbifold twist field 



written as 23 1 



24 



25 



the OPE can be 



dX{z)ag{w,w) '--^ {z — w) '"^ 'Tg{w,w), 

dX{z)ag{w,w) ~ (z — w)^ Tq{w,w), 



where X{w, w) is the intersection point of the two D-branes, (Tg{w, w) is the open string twist field, 



and T and r' are excited twist fields. The local monodromy conditions are then given by 23] 

dX{e^''\z - w)) = e^^'dX{z - w), 

dX{e^^\z - w)) = e-^^'dX{z - w). (3) 

Since we are interested in a four-point interaction, we demand that the four strings are in 
a bounded area such that ^0j = 2, and it is required that the four twist operators ae^{zi,Zi) 
are present at the D-brane intersections. Due to invariance under S'L(2,M), we can set zi = 0, 
Z2 = X, Z3 = 1, and Z4 = Xqq. The field X has a classical piece X^i and a quantum piece Xqu, 
so the interaction amplitude can be factorized into a classical solution resulting from worldsheet 
instantons and a quantum contribution resulting from quantum fiuctuations as 

Z=Y,(^~''''Zgu, (4) 

where 

Sci = ^ f cPz{dx,idXd + dx,idx,i). (5) 

A. The Classical Contribution 

The classical field can be expanded by OPE as 

dXciiz) = a Lo{z), dXci{z) = a J{z), 

dXd{z) = buj'{z), dX,i{z) = bO{z), (6) 

where 

u;iz) = l[{z-x.r^'-'^\ oj'{z) = Yl{z-x,)-'^, (7) 

i i 

a, a, 5, and b are complex constants. On a torus with brane intersection considered [23i], the local 
monodromy condition is given by 

X{e^^'z, e-2-*z) = e^^'^X + (1 - e^^'%U + v), (8) 

where / is the intersection point and v is the lattice translation of the torus. Then from the global 



monodromy conditions we have 23|] 

Ac,X,; = 4e--*(^»-^»+i) sinvr^, sin7r0,+i(/,+i - f^ + v,) = <f {dzdX,i{z) + dzdX,i{z)). (9) 




FIG. 1: Definition of the indices of the vertex and D-branes. 



In Fig. [H the index i denotes the i-ih vertex, /i is the intersection point between stack a and b 
D-branes and so on, and we define Lj,t+i to be the vector along the direction from fi to /j+i- Thus 
vi = qi\Iab\Li2, etc, where qi is an integer and lab is the intersection number. From the closure 
of the quadrilateral ^Vi = 0, we can count the areas from the solutions of the linear diophantine 
equations [23(] 



(10) 



In the case having two independent angles which we will consider later for a specific model, if 
stacks a and c are parallel to each other, then the above diophantine equations can be simplified, 





23l | these Qi can be parameterized 



and the four integer parameters qi are not independent. From 

by two independent integer variables, which we will use in the following area analysis by applying 

a different approach. 

Refer to |23l | for the details of the calculation involving hypergeometric functions from the 
integral of uj and uj', we can write down the minimum classical action when there is non-zero 
worldsheet area on one subtorus as 



23] 






i2 
cl rain 



1 / sin vr^i sin 7r04 v\^ sin 7r02 sin vr^s ^Ij 



+ 



(11) 



27ra' \ sin(7r6'i -F 7r6'4) 2 sin(7r6'2 + vr^a) 2 
where the equation in the parenthesis is exactly the formula for the area of a quadrilateral with 



angles 9i and two sides ^14 and ^723. We will then use this information to calculate the classical 
contribution of the four-point functions. 

B. The Quantum Contribution 



The quantum contribution of the four-point function is given by the correlator 

{ae^ {Zi)ag^ {Z2)crg.^ {Z3)(7g_^ (24)) . 

From the OPE of the stress tensor T{z) with the twist fields we have 



T{z)ag, (zi) 



hi 



+ —. — ^ + 



{Z - ZiY (Z - Zi) 

where hj = \Oj{l — 6j) is the conformal dimension. T{z) also has a relation to X as follows 

- Id.xd^x ~ , ^ ., + T{z) + ■■■ . 

2 [z — wY 

Therefore, we can write 

{T{z)ao^ {zi)ag^ {z2)ag.^ {z3)ag^ (24)) 



lim 



(^ei{zi)(Tg2{z2)ag.^{zs)ag^{z4)) 

{-^dzXd^^Xag^{zi)ae2iz2)crg^{z3)ag^{z4)) 



1 



(crgj {zi)ag2 {z2)(Tg.^ {Z3)ae4^ (za)) 
We may then define the Green's function as 



{z — wY 



g{z,w;zi) 



\dzXdyjXag^{zi)ag2{z2)crg3iz3)ag^{z4)) 



{ag^{zi)ag2{z2)ag^{z3)ag^{z4)) 
Prom Eq. p^ we can find the asymptotic properties oi g{z,w,Zi): 

g{z,w;Zi) ~ -^ + finite for z ^>- w, 

[z — wY 

1 

for z ^ Zj, 

Similarly from Eqs. ([T|) and ([7]) we can write down the equations for X and X: 

dX{z) ~ u'g^z) = \{{z - Xi)-'^ = u;i-eAz). 



{z- 


r, 1 11111 

- wY 
1 


(z- 


-Zi)-^^ 

1 
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Then g{z,w; Zi) can be expanded in the following form 

{Z - Zi){z - Zj)Y{^{w - Zk) 

{w — Zi){'w — Zj){z — wY 



g{z,w;zi) =ujg^{z)ujg^{z)l ^, 



+ A 



(12) 



(13) 



(14) 



(15) 



(16) 



(17) 



(18) 



(19) 



The coefficients a^- can be fixed by the above asymptotic relations, so finahy we can determine 

(r(z)o-e^ (2:1)0-03 (z2)o-e3(^:3)q"e4 (^4)) 
(crei {zi)cr0^ (^2)0-63 {z3)cre4 (^4)) 

= -^E^^%7 ^ T + jE«^.(^- + ^)' + TTr^- (20) 

2^^ ' [z - Zi){z - Zj) 2^ ^^--Zj ■z-'^i ll(-2--2i) 

The quantum part of the monodromy conditions is independent of the contours, so that 



Ac^Xgu = 0= f dzdXgu + f dzdXgu- (21) 

JCt JCi 

We can use this condition to determine A and define two homology cycles Ci and C2 by cutting 
the complex plane between zi and Z2-, and between Z3 and Z4. In this way, we find 

dzg{z,w) + i dzh{z,w) = 0, (22) 

where /i(z, w) is the auxiliary correlation function 

.,- . {-\d^Xdu,XaQ^{z{)(JB.,{z2)(JQ^[,z^)(JB^{z^)) 

h{z,w;zi) = -— = BuJi-eAz)^e,{z)- (23) 

(cre^ (zi)o-e2 (z2)o-03 (23)0-94 (Z4)) 

After using the SL{2,R) invariance, we obtain that Eq. (I2ip turns out to be 

B (p io' {z)dz + A (p uj{z)dz = XoD f y_^aj4(z — 2;j)a;(z)(iz. (24) 

J Gj •y i^i J i^i ^ 

Solving A, the correlator {(Tg-^(Tg^aQ,j^aQ^) then has a form [23|] 



where I(x) is a function of x, 9i, gamma functions of Oi, and hypergeometric functions 23 1. 

Again, we are more interested in the case with two independent angles, therefore by setting 
9i = 1 — 62 = V and ^4 = 1 — ^3 = A, the quantum part of the four-point function can be written 
as 



22 



23| 



Z4g = 167rtx-''(^-'')(l - x)-'^^/(x)-|, (26) 

where 

l{x) = (1 - x)^^-^-^^ {B{y, A) 2Fx{v, A, z^ + A; 1 - x) 2-^1(1 - i^, 1 - A, 1; x) 

+B(l-i/,l-A) 2Fx{\-v,\-\2-v-\;\-x) 2Fi(z^,A,l;x)], (27) 

and 2F\ is the hypergeometric function; B{y^ A) = V {y^Y {X) jV {y + A). We have used the fact that 
2Fi(z^, A, l;x) = (l-x)(^-'"~^) 2F\(\-v,\-\\;x). Note that the quantum contribution Eq. ([26]) 
remains invariant under the interchange of the interior and exterior angles i^, A <-> 1 — z^, 1 — A. 



The general form for a four-point amplitude on T^ with two independent angles is then given by 



Za = 167r 



5/2 



n 



-Vj{l-Vj) 



A 



(l-x,r^^^/,(x,)-V^^exp(-3 



2i:a' 



(28) 



j k,l 

By taking the limit of x, Z^g is fixed to a constant and we obtain the three-point amplitude for 



the Yukawa couplings 
Z3 



m 



M 



2,n 



167r2r(l - Uj)T{l - Xj)T{uj + Aj) 1 



r(i.,)r(A,)r(i 



A.) 



^^M 



A. 



2TTa' 



(29) 



As we will later encounter in the model considered in the next section, where only the D-branes 
on the second torus form closed areas either for four-point or three-point amplitudes, the quantum 
contributions from other tori are just constants, which are then able to be absorbed into the VEVs 
of the Higgs fields. Using a computer code, we find that the ratio of Z^g to Z^g is around 0(10) 
which can then also be absorbed into the VEVs for simplicity. Therefore, we can merely focus on 
the classical contributions to the SM fermion masses and mixings. 

III. A WORKING EXAMPLE OF THE SM FERMION MASSES AND MIXINGS 

A. The Pati-Salam Models 



Let us first review the Pati-Salam model discussed in [17[, where its D6-brane configurations 
and intersecting numbers are presented in Table [H To explain the electron and muon masses in 
this model, we are looking for four-point interactions such as 



i'ab(l^ia(l>b'c(Pbb' or 



'>lb(l^ca(Pcc'(Pbc' > 



(30) 



where (/)Jj^ are the chiral superfields at the intersections between stack a and /? D6-branes. 



stk 


N 


(ni, li){n2, /2)(n3, I3) 


a 


8 


( 0,-1) (!,!)( 1, 1) 


b 


4 


( 3, 1) ( 1, 0) ( 1,-1) 


c 


4 


( 3,-1) ( 0, 1) ( 1,-1) 



TABLE L D6-brane configurations and intersecting numbers for the model presented in [ITj, where the SM 
fermions and Higgs fields are from the intersections on the first torus. This model is constructed from Type 
IIA TV(Z2 X Z2) orientifold. 



From Table HI considering the wrapping numbers on the first two-torus, we find that the 0,R 
image b' of stack b is parallel to stack c, and stack b is parallel to the Oi? image c' of stack c. It is 



10 



required that b and c' or c and b' must overlap to form a closed area, which is not exactly forming a 
quadrilateral. After carefully studying the possible four-point functions on the first torus, we find 
that we can obtain the correct electron and muon masses. However, the techniques for calculating 
such kind of geometric structure for four-point functions have not been clarified yet. Therefore, we 
turn our attention to the other models. To avoid the generic rank one problem, we require that 
the three SM fermion families arise from the intersections on the same two-torus. The number 
of the Higgs bidoublets from that torus is a multiple of 3. In addition, we require that the QR 
image of the U{2)l D-brane stack forms a closed quadrilateral with the f/(4), U{2)l, and f7(2)/j 
stacks to get rid of exotic particles. And there may also arise additional constraints on the D-brane 
wrapping numbers to form a bounded area. Generally a three-generation model without exotic SM 
type particles has a structure like that of a stack of D-branes which lies on the orbifold while the 
other two are images of each other and the ratio of the wrapping numbers is three, for example 
{n,m) = (1,3), just as the model discussed in [171]. However, similar to the above model, we find 
that such types of models cannot form a normal bounded four-point area. Thus, we will consider 
the Model TI-U-3 in [16]. We present its D6-brane configurations and intersecting numbers of the 
observable sector in Table [III &iid the particle spectrum of the observable sector in Table [Till where 
S^j^ and S]^ are the SM singlets from the intersections of stack b and its ili? image b' . 



stk 


N 


(m, li)(n2, l2){n3, h) 


A 


S 


b 


b' 


c 


1 

c 


d 


06 


a 


4 


( 1, 1) ( 1,-3) ( 1, 0) 








3 


0(1) 


-3 


0(3) 


-3 


0(3) 


h 


2 


( 2, 0) ( 1, 3) ( 1,-1) 








- 


- 


6 


0(3) 


6 


0(3) 


c 


2 


( 1,-1) ( 2, 0) ( 1, 1) 








- 


- 


- 


- 


0(1) 


0(1) 



TABLE II: D6-brane configurations and intersecting numbers in the Model TI-U-3 in '16'], where the SM 
ferniions and three of the six Higgs fields are from the intersections on the second two-torus. This model is 
constructed in the supersymmetric AdS vacuum on Type IIA T^ orientifold with flux compactiflcations. 

The superpotential which includes the trilinear Yukawa couplings is given by 



(31) 



where Y^^j^, Yf-^^, and YK are Yukawa couplings, and Q*, U^"^ , -D^*, L* and E'^^ are the left-handed 
quark doublet, right-handed up-type quarks, right-handed down-type quarks, left-handed lepton 
doublet, and right-handed leptons, respectively. The superpotential including the four-point inter- 
actions is 



W4 ~ ^ {Yl-,iQW^H'^Si + Ylf,,Q'D^^H'^Si + Yl],p E^^ H'^ S^) , 



(32) 
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Quantum Number 


Qa 


Q2L 


Q2R 


Field 


ab 


3 X (4, 2, 1,1) 


1 


-1 





Fl{Ql,Ll) 


ac 


3 X (4, 1,2,1) 


-1 







Fr{Qr,Lr) 


1 
ac 


3 X (4, 1,2,1) 


1 







$» 




3 X (4, 1,2,1) 


-1 





-1 


^^ 


fee 


6 X (1,2,2,1) 





1 


-1 


HI, Ha 




6 X (1,2,2,1) 





-1 






fe'c 


3 X (1,2,2,1) 





-1 


-1 


H'u, H'a 




3 X (1,2, 2,1) 





1 






fe'fe 


6 X (1,1,1,1) 





2 





SI 




6x(l,T,l,l) 





-2 





S'l 



TABLE III: The chiral and vector-like supcrficlds in the observable sector, and their quantunr numbers 
under the gauge symmetry SU{A)c x SU{2)l x SU{2)r. 

where Y-%,, Y-'f,,, and Y-'L, are Yukawa couplings of the four-point functions, and Ms is the string 



ijkl'> ijkl^ 



" ijkl 



scale. 



B. Yukawa Couplings from Three-Point Functions 

In order to calculate the trilinear Yukawa couplings, we first note the intersection numbers on 
three two-tori 



l'^ = -1 


r(2) _ o 
-'afe - -^ ' 


r(3) _ _i 

-'afe - -^ ' 


/£) = 1 , 


-'ea ^ ' 


/(3) _ _l 


4^^ = -1 


r(2) _ o 
-'fee - "^ 


/(3) - 2 

-'fee ^ • 



(33) 
From these, we find that d^'^' = g.c.d.{I^^^ , ^J , lea ) = 3, d^^' = 1, and d^^' = 1. The parameters 

(i) 

of the theta functions in terms of the intersection numbers, the brane shifts e^ , the Wilson line 
phases 6a , and the Kahler moduli J*-*' are 11 1 



5(2) = !i!l _ ^ _ ^ + ei^) - 6^ - 6i^^ £(2) 

3 3 3 3 3 

US) _f(3) _ J3) , 2^(3) 



,(3) 



(34) 
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^(3) = _^(3)_^(3)+2ep)^ 



(35) 



K 



(1) 



J(l) 



K 



(2) 



3J(3) 



K 



(3) 



2J(3) 



oc a: 

For convenience we redefine the shift on each torus as 



a' 



(36) 



eW^ 



.1^) + 6, 



(1) 

6 



J2) 

g(i) e(2) = fi_ 



,(2) 



g(2) 
fca 



g(3) = 



tc 



6?) + 2ei=^) 



(37) 



then the theta function of each torus can be written as [111] 



^ 



(/)(-) 



,(r)^ 



7ri(<5('-)+/r)2K('') 27ri(<5('-)+/r)</.('-) 



K^-^) = J2e'--'-^ -'-^ e 



(38) 



/rSZ 



where r = 1,2, 3, so the Yukawa couphng constant can be expressed as 



Y{ijk} = Zsgaabc I _[ ^ 



r=l 



Sir) 



(kM) 



(39) 



where Z^q stands for the quantum contribution to the instanton amplitude, and aabc = 

n.sign(/l;)/£)/(;;)). 

For the first torus the intersection number is one, so its contribution is just a constant. On the 
second torus v"^' ,j^'^> ,k^'^' range from — 2, and on the third torus k^"^' goes from to 1. Thus, 
although the total number of Higgs state is k^"^' x k^^' = 3x2 = 6, only three linear combinations of 
six Higgs fields can provide the SM fermion Yukawa couplings since the intersection number between 
h and c stacks of D6-brane is 3 on the second two-torus. Explicitly, there is only one SM Higgs 
linear combination from the two bidoublet fields {k^"^' ,k^'^') = {(A;*^^\0), (A;'^', 1)}. Furthermore, 
since the triplet of intersections is connected by an instanton, the selection rule for the indices of 
the second torus 



■(2) ^ -{2) ^ ^(2) ^ Q ^^^ 3 ^ 

should be satisfied. Then we can choose the Yukawa coupling matrices as the following form 

/^00\ /00S\ /oC7o\ 

A 



(40) 






C7 

\ s y 






' ^k=2 



\c y 



BOO 

\ ^y 



(41) 
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where 



A = ^ 



.(2) 



0(2) 



a' 



.(2) 



+ 



</.(2) 



(^), c. 



a' 



e(2)_ 
,/.(2) 



,3J(2), 



a' 



(42) 



where we simply set here s'^^ = k^"^'. Similarly in the third torus, there is only one parameter 



y(3) ^ 



i9 



e(3) 

0(3) 



.2J(3) 



a' 



A,, y£\~^ 



=(3) 



+ 



0(3) 



2/(3) 



a' 



(43) 



Therefore the classical part of this three-point couplings is given by 



/ 



Zscll 



AA3 ^ 

CAs 
BA3 



A 



Z3d2 



Z'iclA 



\ 

I ADz ^ 

CDs 

V BD3 y 



■^3ci3 



-^3c/5 



/ 



^3c«6 



/ 


CA3 


\ 


BAi 








\ 





AA3/ 


/ 


CD-i 


\ 


BD-i 








\ 





ADi) 



(44) 



'^0 BAi 

AA3 
VCA3 ^ 

'^0 BD^ 

^L»3 
yCDs 

If each of the Higgs states resulting from intersections localized between stacks h and c develops a 
VEV, vf, where z = 1, • • • ,6 for each of the six Higgs states, and <I> is an index for (n)p-type quarks, 
((i)own-type quarks, and (/)eptons, respectively, the total effects of Yukawa couplings including the 
quantum part will be 

U{A;vf + D^vf) C{A^vt + D^vt) B{A;vt + D;vf)^ 
Zi = Zs, BiA^svf + D^vi)A{Asvf + D^vf)CiA^svf + D:,vf) , (45) 

\C{A3vt + Dsvf) BiAsvf + Dsvf) A{A,vi + D,vf) j 

where v^ = v\ . Thus, it is clear that only three linear combinations of the six Higgs states 
contribute to the Yukawa couplings: A^Hf + D^Hf, A^Hf + D^Hf, and A^Hf + D^H^, where 
for simplicity we neglect the normalization. 



C. Yukawa Couplings from Four- Point Functions 



The formula for the area of a quadrilateral in terms of its angles and two sides and the solutions 
of diophantine equations for estimating the multiple areas of the quadrilaterals from non-unit 
intersection numbers are given in 23] ■ In the present discussion, we are considering a model 



which possesses only two independent angles, so we need at least two parameters to describe all 
quadrilaterals. In addition to these formulae, there is a more intuitive way to calculate the area for 
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these four-sided polygons with only two independent angles. A quadrilateral with two independent 
angles is a trapezoid, and a trapezoid can be always taken as the difference between two similar 
triangles. Therefore, since we know the classical part is 



it is equivalent to write 



Zid ~ e" 



^quad 



■'Ad 



-\Atri—A^ri\ 



(46) 



(47) 




FIG. 2: A picture of two trapezoids with different field orders. The red brane repeats in a next cycle and can 
still form a similar trapezoid with the blue brane. This coupling also contributes to the four-point function. 



Taking the absolute value of the difference reveals that there are two cases: Atri > A[^- and 
Atri < A'^^-, as shown in Fig. [2j From the figure we can see the two trapezoids are similar with 
different sizes, but the orders of the fields corresponding to the angles are different, which is under 
an interchange ol 9 ^^ 1 — 9 , 6 = u, X. These different field orders may cause different values for 
their quantum contributions. However, we have shown above that this angle transformation will 
not affect the quantum contribution, so these two cases are on equal foot to sum up. Therefore, we 
are able to employ the same techniques which have developed for calculating the trilinear Yukawa 
couplings. 

For a trapezoid formed by the stacks a, b, b', c, we can calculate it as the difference between 
two triangles formed by stacks a, b, c and b' , b, c. In other words, they share the same intersection 
Ifjc- Therefore, if we use this method to calculate the trapezoidal area, we should keep in mind 
that the intersection index k for Ibc remains the same for a certain class of trapezoids when varying 
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FIG. 3: A diagram showing the areas bounded by stacks of D-branes which give rise to the Yukawa couphngs 
for quarks and leptons via world-sheet instantons. The Yukawa couphngs of the up- type quarks are from 
the areas by stack ai, 6, ci, the down- type quarks by stack ai, b, C2, and the leptons by 02, b, C2. The 
four-point function corrections to the Yukawa couplings of the up-type quarks are from the areas by stack 
ai, 6, 6', ci, the down-type quarks by stack oi, 5, b', C2, and the leptons by 02, 6, b' , C2- 



other intersecting indices. Here we set indices i for lab, j for lea, 1 for h'b, and j for Icb', as shown 
in Fi..ffl We .nay clcnla.e .he »eas of the triangles a. we did in the t.iiinear Yni^wa couplings 
above [11 1 



A^ 



ijk 



-{2TrfAT,2\IabIbcIca\ {^ + ^ + — + e + l)\ 

^ ^ nh J rn. J hr 



A. 



ijk 



^ J^b'b J^cb' J^bc 



(48) 



where i, j, k and i, j, k are using the same selection rules as Eq. (j40p . Thus, the classical 
contribution of the four-point functions is given by 



^4d 



E 

11 



g 2-it' 'i'' ^ijk 



(49) 



Note that this formula will diverge when Aijj, = A^i^, which is due to over-counting the zero 
area when the corresponding parameters in Eq. (j48p are the same. In such a case, Z^^i = 1 + 
^^/^^^^"z^'*^*" '•''='. We will not meet this special situation in our following discussion. 
In the model of Table HH in addition to the intersection numbers in Eq. ()33p . we have 



^bb' - ^ ■> 


r(2) _ 
^bb' — "* ' 


r(3) _ 2 
^bb' - ^ 


/(^) -n 


j(2) __3 


t{3) _2 




r(2) _ 


j(3) _ 



(50) 
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and the number of the singlets from the SU{2)i anti-symmetric representation is given by M.Anti = 
Wi^-^w "'"-^boe)/^' ^^ ^^ obvious that we have six vector-Uke anti-symmetric fields for SU{2)l, three 
from the second torus and two from the third torus. The matrix elements on the second torus from 



the four-point functions can be written in terms of Ojo-jj as 



/ 



ooooo 







ai200 
12100 

aoiio 
oioio 

\ 02210 J 

00220 

aii2o 

\ 02020 



/ 



ooioi 

oiooi 

\ 02201 J 

' 00211 

oiiii 

\ 02011 

' 00021 

01221 

\ 02121 



/ 



O0202 

O1102 

\ 02002 y 

/ O0012 \ 

01212 

\ 02112 j 

I O0122 » 

O1022 

\ ^ 02222 / 



The relations between these elements are shown in Table IIVI 



Parameters 


5 


d 


fc = 


fc = l 


k^2 


A 








aoooo 


lull 


12222 


B 


1 

3 





12100 


10211 


11022 


C 


-1 
3 





11200 


12011 


10122 


V 





1 

3 


10021 


11102 


12210 


£ 


1 

3 


1 
3 


12121 


10202 


iioio 


T 


-1 
3 


1 
3 


11221 


12002 


lOllO 


Q 





-1 
3 


10012 


11120 


12201 


n 


1 
3 


-1 
3 


12112 


10220 


llOOl 


1 


-1 
3 


-1 
3 


11212 


12020 


lOlOl 



(51) 



TABLE IV: The matrix elements in terms of parameters. 6 



-L+feandd=7^ + y2- 



Note that since we need to sum over I and i, an integer shift has no effect, so ^ is equivalent to 
equivalent to | for d and 6, etc. 



IS 



In this model, we have six SM singlet fields S]^ and three Higgs-like states H'^^^. Similar to 
the Higgs fields H^ ^, only three linear combinations of the six S^ can contribute to the four-point 
Yukawa couplings. Thus, the four-point Yukawa couplings involve three SM singlet fields S)^ and 
three Higgs-like states H'^^. If their VEVs are denoted as Uj and Wj respectively, then the complete 
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contribution will be 

' AuiWi + QU2W3 + 'DU3W2 TuiW2 + TU2WI + CU3W3 SuiWs + BU2W2 + TLu^wi 

ZacI ~ HuiW2 + £u2Wi + Busws Vuiwz + AU2W2 + Qu-iWi Cuiwi + IU2W3 + TU3W2 ■ (52) 

^TuiW^ + CU2W2 + lu^Wi BuiWi + HU2WS + £U2,W2 GuiW2 + VU2W1 + Au'j,Wij 

IV. A NUMERICAL EXAMPLE 
A. The SM Fermion Masses and Mixings at the GUT Scale 

The main reason for the addition of these four-point corrections is to better match the Yukawa 
coupling matrices to those obtained by running the renormalization group equations (RGEs) up 
to the GUT scale. In particular, we would like to better match the lepton masses. It should be 
remembered that the parameters from the theta function and the parameters in Table IIVI which 
depend on the D-brane shift parameters are not independent. Thus, the different elements of the 
Yukawa mass matrix are related, as discussed in [XT\ which strongly constrains the form that they 
may take. As before, we will perform a transformation on the SM fermion mass matrices which 
are obtained from the RGE running of the experimental values for the SM fermion masses up to 
the GUT scale to make the comparison with the theoretical results. If we define D^ and D^ as the 
mass diagonal matrices of the up- and down-type quarks respectively, the transformations are 

UlMJjf = D^, UlMaU^^ = Da, VcKM = UlUt\ (53) 

and the squared mass matrices H^ and Hd are 

Hu = MuMl Hd = MdMl (54) 

For simplicity, we assume that the quark mass matrices M„ and M^ are Hermitian. If we take a 
case in which M^ is very close to the diagonal matrix, or in other words C/^ and Uf^ are very close 
to the unit matrix with very small off-diagonal terms, we have 

VcKM ~ U'^U''^ ~ f/", (55) 

where we have transformed away the right-handed effects and make them the same as the left- 
handed ones. Then the mass matrix of the up-type quarks turns out as 

^u ~ vl^j^j^jDuVcKM- (56) 
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VcKM 



At the GUT scale, the CKM quark mixing matrix is given by [27] 

^ 0.9754 0.2205 

-0.2203e°-°03°i 0.9749 

\ 0.0075e-i9°* -0.0311e^-0°^ 
and Du and D^ are 

/ 0.0000139 o\ / 0.00141 



-0.0026i ^ 
0.0318 
0.9995 y 



(57) 



A, 



mt 



v 








0.00404 
1 



, Dd = mb 



\ 








o^ 

0.0280 
ly 



(58) 



so that the SM fermion mass matrices (we exclude the phases in this discussion) that we would 
like to match are given by 

/ on09fifi 00109 00747 \ /nnni/Li n n\ 



\Mu\ = rrit 



0.000266 0.00109 0.00747 

0.00109 0.00481 0.0310 

\ 0.00747 0.0310 0.999 ) 



\Md\ = Dd = nib 



\ 



0.00141 






0.0280 
1 



(59) 



We use the relation ^^ 



1.58, as well as require the eigenmasses of the leptons to be 
{me, mi,,mr} = m^{0.000217, 0.0458, 1}. 



(60) 



B. The Numerical Results 



As mentioned above, there are six Higgs states in this model, and only three of them can provide 
the SM fermion Yukawa couplings from three-point functions on the second torus. And from Eq. 
()45|) we can see explicitly only three linear combinations of the VEVs dominate the SM fermion 
mass matrices. To fit the diagonal terms in the up-type quark mass matrix, we have a limited 
number of degrees of freedom to work with on the off-diagonal mixings. Thus, it is hoped that the 
four-point corrections will improve this situation. In short, the situation for the Yukawa couplings 



17| due to the highly dependent relation 



of this model is not as ideal as the model analyzed in 
between the diagonal and off-diagonal terms which results from fewer Higgs fields. Therefore, for 
the best fits to the SM fermion masses and mixings at the GUT scale, we will form the Yukawa mass 
matrices such that the diagonal terms result from both the three-point and four-point couplings, 
while the off-diagonal elements arise strictly from the four-point corrections. Since we also expect 
the down-type quark mass matrix to be diagonal, it is necessary to choose reasonably small off- 
diagonal terms for the mixing of the first two quark generations, so we set the Kahler parameter 
30.0. 



3J(2) 
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1. Quantum Contributions 

As mentioned above, the quantum contributions from the three-point and the four-point inter- 
actions can be ignored since they are able to be absorbed into the VEVs of the vector-hke fields: 



Zsgf ; U[w'j = Z^gUiWj 



(61) 



However it is interesting to study how much they may affect the VEVs. The angles between the 
branes on the second torus are 



Avr = Zab = Zbb' = 7r/2 ; i/tt = Zac = Zb'c = tt/A 



(62) 



so from Eqs. (j29p and (j26p we have 



Z^?~ 12.95; zi?(x)U^o.5~120. 



"3g 



^4g 



(63) 



In this case, the plot of Z4q as a function of x in the range x = (0, 1) is shown in Figure [U 
where we see that ^4^ approaches its three-point function limits of the corresponding fields when 
X approaches or 1 as a constant, so we confine our interest of ^4^ near x = 0.5. We find that 

Z4g/Z3g ~ 0(10). 




0.2 



0.4 



0.6 



0.8 



FIG. 4: Ziq as a function of x in the range x = (0, 1). 



2. Yukawa Couplings from Three-Point Functions 



Let us first consider the SM fermion Yukawa couplings only from the three-point functions. For 
the best fits, we set appropriate values of Afv + Dfv to exactly match the diagonal terms of 
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<l>\Parameters 


e(2' (shift) 


A* 


B* 


C* 


Al 


Df 


Alvf+Dfvf 


Alv'i + Divi 


Alvi + Divi 


<1> = u up quarks 





1 


0.000028 


0.000028 


1 





0.000266mt 


0.00481mt 


l.Omt 


$ = d down quarks 


0.0775 


0.567749 





0.002094 


1 





0.00248mf, 


0.0493m6 


1.761m6 


$ = I leptons 





1 


0.000028 


0.000028 


0.901 


~ 


0.00224mf, 


0.0444m6 


1.587m6 



TABLE V; The matrix elements and the parameters when only trilinear Yukawa couplings considered. Note 
that vf — vl. 



the SM fermion mass matrices at the GUT scale in Eqs. ()59p and ()60l) . The parameters and the 
VEVs in the quark and lepton mass matrices are hsted in Table |Vl 

The SM fermion mass matrices from the three-point function contributions are 



/ 



\M- 



3m I 



mt 



0.000266 0.000028 1.34 • 10 



-7\ 



0.000028 0.00481 



1.34 • 10" 



-7 



~0 
1.0 



/ 



\M- 



3d\ 



rub 



0.00141 0.00369 







\ 



0.0280 5.19 • 10-*^ 

0.000103 1.0 



/o. 



IM3/I ~ TTlb 



00224 0.0000444 1.24 • 10 



-6\ 



(64) 



0.0000444 0.0444 ~ 

\^1.24-10"^ ~0 1.568 J 

Therefore, with only three-point functions, we can obtain the correct SM quark masses, and tau 
lepton mass. However, the quark CKM mixings are too small, and the electron and muon masses 
are far from the desired values. 



Yukawa Couplings from Three-Point and Four-Point Functions 



Let us study the SM fermion masses and mixings from both three-point and four-point functions. 
To be concrete, we present the SM fermion mass matrices from three-point and four-point functions 
separately in our best fits. First, we present the parameters and the VEVs for three quark and 
lepton mass matrices of the Yukawa couplings from the three-point functions in Table IVII 



$\Parameters 


e(2' (shift) 


A" 


B* 


C* 


Al 


Di 


Alvf + Dfv2 


Alv$ + Divi 


Alvi + Divi 


^ — u np quarks 





1 


0.000028 


0.000028 


1 





0.0000250mt 


0.000262mt 


0.981mt 


$ = rf down quarks 


0.0775 


0.567749 





0.002094 


1 





0.002820m6 


0.08784mt, 


1.754m6 


<& — I leptons 





1 


0.000028 


0.000028 


0.901 


~0 


0.002541m6 


0.07916mt, 


1.581m6 



TABLE VI: The matrix elements in terms of parameters. Note vf — w'. 
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The SM fermion mass matrices from the three-point interactions are 

/0.OOOO25O 0.0000275 ~ \ / 0.00160 0.00367 



\M: 



3pt| 



3« 



mt 



\ 



0000250 0.0000275 ~0 
0.0000275 0.000262 ~0 
~ <- 0.981 



[M|fl 



■mb 



\ 
0.0499 0.00000591 

\0.000184 0.996 J 



I 



\m: 



3pi| 



3/ 



nib 



0.00254 0.0000443 2.22 • 10 



-6\ 



(65) 



0.0000443 0.0792 7.12 • 10"' 
\^2.22 • 10"*^ 7.12 ■ 10-s 1.581 j 

Second, for the Yukawa coupHngs from the four-point functions, we present the parameters and 
the shifts of the D-branes in Table EUl and the VEVs of three SI and H'^^ in Table Eml 



$\Parameters 


,(2) 


4?' 


A" 


B* 


C* 


P* 


£* 


^ 


e* 


7i* 


J* 


$ = u up quarks 





0.32 


0.000064 


0.679089 


0.679089 


0.000012 


0.670140 


0.670140 


1.335992 


0.000029 


0.000029 


$ = d down quarks 


0.08 


0.03 


0.595678 





0.00257000 


0.000007 


0.028872 


0.001674 


0.000313 


0.000595 


0.073686 


$ = I leptons 





0.03 


0.932776 


0.000031 


0.000031 


0.000004 


0.165161 


0.165161 


0.000171 


0.187111 


0.187111 



TABLE VII: The matrix elements in terms of parameters. 



$\Parameters 


ui{Ms) 


U2(Ms) 


us{Ms) 


wi(mi,) 


vj2{mb) 


W'i{mb) 


^ — u up quarks 


0.08 


0.06 


0.23 


0.0148 


0.18 


0.003 


$ = d down quarks 


0.08 


0.06 


0.23 


-0.004 


-0.612 


0.03 


$ — I leptons 


0.08 


0.06 


0.23 


-0.004 


-0.612 


0.03 



TABLE VIII: The VEVs of three Si and iJ^\ 



The SM fermion mass matrices from the four-point interactions are 



/ 



\mZ\ 



mt 



0.000241 0.00106 0.00750 



\m: 



4pt< 



il 



rrib 



0.00106 0.00455 0.0285 
\ 0.00750 0.0285 0.0192 J 

/-0.000299 -0.00920 0.000223^ 

-0.00920 -0.0343 -0.0229 
\ 0.00657 -0.0229 0.00623 y 



<'l 



nib 



/^ -0.000191 -0.00359 0.0000688^ 

-0.0000361 -0.0219 0.000104 
\ 0.000158 -0.00406 0.00409 J 



(66) 
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Therefore, we sum over the two contributions, and obtain 



/o. 



IM,, 



mt 



000266 0.00109 0.00750 



/ 0. 



\Mi\ ^mb 



0.00109 0.00481 0.0285 
\ 0.00750 0.0285 1.0 J 

( 0.00224 -0.00915 0.000223^ 
-0.00915 -0.0449 -0.0229 
0.000223 -0.0229 1.587 



|Md! 



rah 



00141 0.0000828 0.0000688 



\ 



-0.0000361 0.028 -0.0000979 
\ 0.0000258 -0.00406 1.0 j 



\ 



(67) 



and the mass eigenvalues of the leptons are 



{me,m^, rrir) = m5{0.000348, 0.0464, 1.588} = m^{0.000219, 0.0292, 1}. (68) 

Comparing with the SM fermion masses and mixings in Eqs. (j59p and (j60p , we find that the off- 
diagonal terms of the up-type quark mass matrix are basically from the four-point contributions, 
and it is quite difficult to fit exactly the same off-diagonal terms due to the four-point function 
matrix structure. For the down-type quark matrix, we still cannot eliminate the off-diagonal terms. 
However, we can suppress them such that they are in an acceptable range. In addition, the quark 
CKM mixing matrix is 

/ n Q77 919 nnmnQx 



VCKM 



0.977 0.212 0.00109 



(69) 



0.212 0.977 0.0298 
\0.00738 0.0289 0.9996/ 

and then there are some deviations for the quark mixing terms in Eq. (j57p . but they are in an 
acceptable range. Finally, we are able to decrease the electron mass eigenvalue to the correct value 
by the off-diagonal terms from the four-point contributions, but the muon mass eigenvalue is still 
about 36% smaller than the desired value. 



CONCLUSION 



We have discussed corrections to the SM fermion Yukawa couplings in intersecting D6-brane 
models due to four-point interactions and presented a working example, and demonstrated that 
these corrections can improve the best fits for the SM fermion masses and mixings. In a concrete 
model, we first calculated the SM fermion masses and mixings from three-point functions. Consid- 
ering only these contributions, we can obtain the correct quark masses and tau lepton mass, but 
the CKM quark mixings are not large enough, and the electron and muon masses are far from the 
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desired values. After including the corrections to the SM fermion Yukawa couplings from four-point 
functions, we can obtain the correct quark masses and CKM mixings, and the correct electron and 
tau lepton mass scales. However, the muon mass is still around 36% smaller than the desired value. 
In this work, we have not considered the moduli stabilization problem, and have essentially 
treated the moduli VEVs as free parameters. As is obvious, there exists fine-tuning in our discussion 
of the D-brane positions and the Higgs VEVs. The four-point interactions have nine dependent 
worldsheet instanton parameters as well as six additional string-scale Higgs particles and their 
VEVs. To stabilize these undetermined variables, one may consider a Z2 x Z2 orientifold model 
with discrete torsion where the D-branes wrap rigid cycles, thus stabilizing the open-string moduli. 
In this case, one would also expect corrections to the Yukawa couplings from D-brane instantons. 
Indeed, E2-branes required for this construction must also wrap rigid cycles, so there is additional 
motivation to consider this background. We plan to pursue these possibilities in our future research. 
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